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(Dated: August 8, 2008) 

We present a measurement of the cross section for Z boson production times the branching fraction 
to tau lepton pairs o{pp — > Z + X)-Br(Z — > t + t~) in pp collisions at y/s = 1.96 TeV. The measure- 
ment is performed in the channel in which one tau lepton decays into a muon and neutrinos, and the 
other tau lepton decays hadronically or into an electron and neutrinos. The data sample corresponds 
to an integrated luminosity of 1.0 fb _1 collected with the DO detector at the Fermilab Tevatron Col- 
lider. The sample contains 1511 candidate events with an estimated 20% background from jets or 
muons misidentified as tau leptons. We obtain a ■ Br = 240 ± 8 (stat) ±12 (sys) ±15 (lum) pb, which 
is consistent with the standard model prediction. 
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PACS numbers: 13.38.Dg, 13.85.Qk 



The resonant production of tau lepton pairs is as inter- 
esting for the study of standard model (SM) physics as 
the production of lighter lepton pairs. For new phenom- 
ena, especially for decays of particles coupled to mass, 
such as SM or supersymmetric Higgs bosons, the de- 
tection of resonant pairs of tau leptons becomes even 
more interesting. This is due to the fact that tau lep- 
tons are much heavier than the other leptons, increasing 
the chance that these new phenomena would be observed 
first in this channel. Unfortunately the detection of tau 
leptons is far more difficult than that of muons or elec- 
trons. 

A measurement of <j(pp — * Z + X)-Bi(Z — > t + t~) in 
pp collisions at y^s = 1.96 TeV is described in this Let- 
ter. The analysis is based on an event sample containing 
a single isolated muon from a tau lepton decay and a tau 
candidate reconstructed as a narrow jet that could be 
produced by a tau lepton decaying either hadronically or 
into an electron and neutrinos. This measurement is of 
interest not only as a test of the SM prediction but also 
because any excess over the expected a-Br could be an 
indication of a source other than Z bosons for events con- 
taining tau lepton pairs, such as the Higgs boson The 
precision of this result is significantly improved compared 
to earlier publications [2], |3j ■ 

The analysis presented here is based on data collected 
between September 2002 and February 2006 by the DO 
experiment, corresponding to an integrated luminosity of 
1003 ± 62 pb" 1 [|. 

The DO detector 5] is a general purpose, axially 
and forward-backward symetric detector, consisting of a 
central-tracking system located within a 2 T supercon- 
ducting solenoidal magnet, surrounded by three liquid- 
argon/uranium calorimeters and a muon detector. The 
spatial coordinates of the DO detector are defined using 
a righthanded Cartesian system with the origin in the 
center of the detector. The positive z-axis is the direc- 
tion of the proton beam, the positive y-axis points up- 
wards and the positive x-axis points out of the Tevatron 
ring. The azimuthal angle </> is measured with respect 
to the positive x direction. Pseudorapidity is defined as 
7] = — ln[tan(#/2)], where the polar angle 9 is measured 
with respect to the positive z direction. The tracking 
system has coverage up to t) w 3. The calorimeter con- 
sists of a central section (CC) covering \i]\ < 1.1 and two 
end calorimeters (EC) that extend coverage to |7/| ~ 4.2, 
all housed in separate cryostats and segmented into cells 
of dimensions 0.1 x 0.1 in?] — 4> space [6|. The muon 
system Q provides a coverage up to rj w 2 and is lo- 
cated outside the calorimeter; it consists of a layer of 
tracking detectors and scintillation trigger counters be- 
fore 1.8 T iron toroids, followed by two similar layers 
after the toroids. Luminosity is measured using plastic 



scintillator arrays located in front of the EC cryostats, 
covering 2.7 < |?7| < 4.4. A three level trigger system is 
designed to select most interesting events based on pre- 
liminary information from the tracking, calorimetry, and 
muon systems, reducing the number of recorded events 
from the collision rate of w 2 MHz to a rate of « 50 Hz, 
which is written to tape. 

The triggering strategy used in this analysis is based 
on the tau lepton which decays into [iv^v-,-. A single 
muon trigger requiring hits in the muon system in com- 
bination with a high transverse momentum (pt) track 
reconstructed in the central tracking system is required. 
The average trigger efficiency ultimately parametrized as 
a function of </3, rj and z using a data sample of Z — > /-< + /i _ 
events, is (52.3 ± 1.4)%. No dependence on the muon pt 
is observed above 15 GeV. 

Most backgrounds as well as the efficiency of the se- 
lection for signal Z — > t + t~ events are estimated using 
Monte Carlo (MC) simulations. All simulated samples 
are generated with PYTHIA [8J] using the CTEQ6.1L par- 
ton distribution function (PDF) set. Simulation of the 
DO detector is done using GEANT3 . Noise in the detec- 
tor and the contributions from other simultaneous inter- 
actions are simulated by adding random untriggered data 
events to the MC simulation. These events were chosen 
such that the effective instantaneous luminosity distribu- 
tion in MC is the same as in data. The code used for the 
reconstruction of simulated events is identical to the one 
used for data. 

Corrections are applied to all MC events to obtain 
overall good agreement between the simulation and col- 
lider data. The momentum scale and resolution for 
muons in the MC are tuned to reproduce the Z boson 
invariant mass distribution observed in data. Similarly, 
the jet energy resolution is tuned to match that observed 
in data for each region of the detector. The pr spectrum 
of the Z boson for events generated with pythia has a 
different shape than that measured in data; therefore the 
Pt of the Z boson is reweighted to fit the direct mea- 
surement in data [Io| . Small differences in acceptance 
between data and simulation are corrected for by weight- 
ing the simulated z position of the primary vertex in MC 
events to reproduce that observed in data. 

Reconstruction efficiencies for muons and tracks are 
calculated both in data and MC using samples of Z — > 
events. Efficiency correction factors for MC events 
as a function of muon or track </3, r\ and z are applied. 
The signal or background samples are normalized to the 
expected number of events evaluated using the luminos- 
ity of the data sample and the theoretical values of the 
next-to-next-to-leading order (NNLO) cross sections in 
the case of Z boson production [ll|, EH or next-to- leading 
order (NLO) cross sections for all other processes where 
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FIG. 1: NN output distributions for (a) type 1, (b) type 2, and (c) type 3 tau candidates. The ratio of signal to background 
is arbitrary, but the relative amounts of type 1, type 2, and type 3 events in background and signal are not. The distributions 
are normalized with respect to each other such that the sum over the three types is 1 for both signal and background. 



the NNLO calculation is not available. W boson produc- 
tion is normalized from data. 

In this analysis, muons are identified starting with 
their signature in the muon detector. The track recon- 
structed from hits in the muon layers is required to match 
a track from the central tracking detectors. The muon 
momentum is measured using only the central tracking 
detectors. 

A tau candidate is a collection of (i) a calorimeter clus- 
ter reconstructed using a simple cone algorithm [l3| , (ii) 
tracks associated with the calorimeter cluster of which 
at least one has px > 1.5 GeV but with a total in- 
variant mass less than 1.8 GeV, and (iii) electromag- 
netic (EM) sub-clusters constructed from the cells in 
the EM section of the calorimeter. The size of the 
cone used for reconstruction of the calorimeter cluster is 
7Z = y/ (A0) 2 + (Ary) 2 = 0.5, where Acj> is the difference 
in azimuthal angle, and Arj the difference in pseudora- 
pidity between the cone axis and each of the calorimeter 
towers. Isolation variables are calculated using a cone of 
1Z = 0.3. The tracks associated with the tau candidate 
must also be contained within this 1Z — 0.3 cone. 

Tau candidates are classified as type 1, 2 or 3, depend- 
ing on the numbers of tracks and EM sub-clusters they 
possess. Type 1 tau candidates have exactly one associ- 
ated track and no EM sub-clusters, type 2 have one asso- 
ciated track and one or more EM sub-clusters, and type 
3 have at least two associated tracks. These categories 
correspond roughly to pure one-prong decays, one-prong 
plus neutral pion decays as well as decay into electrons, 
and three-prong decays of the tau lepton. 

Due to the large number of jets reconstructed as tau 
candidates, additional selection criteria must be applied 
in order to distinguish tau leptons from jets. Three neu- 
ral networks (NN), one for each tau type, are trained 
using Z — > t + t~ MC events as signal and events with a 
jet back-to-back with a non-isolated muon from data as 
background. The NNs use isolation variables based on 



tracks, hadronic and EM calorimeter clusters, as well as 
shower shape variables and correlation variables between 
calorimeter and tracks. Figure [T] shows the discrimina- 
tion obtained using the NNs. Requiring that the NN out- 
put be larger than 0.9 results in a background rejection of 
almost a factor of 50 for all three tau types. This reduces 
the probability for a jet to be misidentified as a tau lep- 
ton to 1.1% for the sum of all types (from 52% without 
the NN output requirement), while maintaining a total 
efficiency of close to 70% for tau leptons which decay 
hadronically or to an electron and neutrinos. Electrons 
are treated as type 2 tau candidates, since the efficiency 
for them to be reconstructed as such and pass the NN 
output requirement is 98%. For a complete description 
of the neural networks and details on their performance 
see Ref. 

The variable chosen to best illustrate the Z — > t + t~ 
signal is the visible mass, given by: 

Visible Mass = ^(P^ + P r + Pr) 2 , (1) 

where P MiT = (E^ T ,p^ T ,p v ^ T ,p^ T ) are the four- 
momentum vectors of the muon and the tau candidate, 
and pr = ($t, $ti $t^)i w i tn $t being the missing 
transverse energy in the event and Tpj, , f T being its pro- 
jections along the x and y directions. The uncorrected 
missing transverse energy is defined as the vector equal 
in length and opposite in direction to the vectorial sum 
of transverse energies of the calorimeter cells. The trans- 
verse momenta of muons are subtracted from this vector, 
after corrections for the energy deposited by the muons 
in the calorimeter have been applied. When the tau can- 
didate matches a reconstructed electron, the energy cor- 
rections derived for electrons are applied. For jets cor- 
responding to tau candidates, the tau energy corrections 
described below are applied. Jet energy corrections ap- 
plied to all other jets in the event are propagated to the 
missing Et calculation. 
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To compare the visible mass distributions of the tau 
pairs between data and MC, it is important to have the 
correct energy scale for the tau candidate. For type 1 
tau candidates, the momentum of the track is used as the 
best estimate of the energy of the tau candidate when the 
tracking resolution is superior to the calorimeter energy 
resolution (up to calorimeter cluster energy of 70 GeV). 
For type 2 candidates matching electrons, the energy cor- 
rections derived for electrons are applied. For type 2 can- 
didates not matching electrons and type 3 tau candidates, 
the energy is estimated using 



E c 



,Pi 



trk 



trk 



(2) 



where p* rk is the momentum of track i associated with 
the tau candidate, E c&[ is the energy deposited by the 
tau candidate in the calorimeter, and i?(p* rk , rj) repre- 
sents the response of the calorimeter to the 7r ± which 
produced track i associated with the tau candidate, as a 
function of the energy and rapidity of the . Typically, 
0.6 < R(pl rk , rj) < 0.9. As the resolution of the calorime- 
ter is better than that of the tracking at calorimeter 
cluster energies higher than 70 GeV (type 1), 100 GeV 
(type 2), or 120 GeV (type 3), the energy of the calorime- 
ter cluster is used in these cases, after applying rj and 
energy dependent corrections obtained from MC. 

The default program in the DO GEANT simulation 
for hadronic interactions, geisha 15j, does not repro- 
duce the charged pion response in data well. Therefore 
gCALOR is used for a more precise simulation of sin- 
gle charged pion interactions. The charged pion response 
obtained using these special simulations was found to 
be in reasonable agreement with preliminary data mea- 
surements in the central calorimeter [Tt| . The energy 
measurement for neutral particles, mostly important for 
type 2 taus, is dominated by electromagnetic showers in 
the calorimeter. The simulation of electromagnetic show- 
ers in GEANT is sufficiently accurate for the purpose of 
this measurement. 

The preselection requires one isolated muon recon- 
structed within the pseudorapidity interval \r)\ < 1.6. 
The transverse momentum of the muon as measured 
by the central tracking detectors must satisfy p T > 
15 GeV. No other muon matched to a central track with 
Pt > 10 GeV is allowed in the event. The muon iso- 
lation requires the sum of energies of all cells situated 
in a hollow cone around the direction of the muon with 
0.1 < 1Z < 0.4, as well as the sum of all tracks in a 
cone of 1Z < 0.5, excluding the muon track, to be less 
than 2.5 GeV. 

The preselection further requires one tau candidate 
with pt > 15 GeV, \r)\ < 2, scalar sum of the transverse 
momenta of all tracks associated with the tau candidate 
> 15 GeV for types 1 and 3 and > 5 GeV for type 2 tau 
candidates, NN > 0.3, and no other muon matching the 
tau candidate. Type 3 tau candidates with two tracks 



are only considered if both tracks have the same charge. 
The tau candidate is required to have a charge with op- 
posite sign to that of the muon. The distances in the z 
direction at the track's point of closest approach between 
the muon and the primary vertex, the tau candidate and 
the primary vertex, as well as the distance between the 
muon and the tau candidate must be less than 1 cm. 

In total 8316 events pass these criteria. To reduce the 
W + jets and the Z — > backgrounds, another se- 

lection criterion is used, based on a variable which gives 
an approximation of the W boson mass, referred to as 



2£ I ,£ A1 (l-cosA</>), 



(3) 



where E v = f!j> • E^/p T is an approximation of the neu- 
trino energy, and A<j> is the angle between $t and the 
muon in the transverse plane. 

For the final selection, all the preselection criteria are 
applied. Additionally, the lower limit on the NN output 
for the tau candidates is raised to 0.9 for types 1 and 2, 
and to 0.95 for type 3 tau candidates. The final selection 
also requires m* < 20 GeV. A total of 1511 events pass 
all the selection criteria. 

The dominant remaining background arises from mul- 
tijet processes, mainly from bb events where the muon 
isolation requirement is met and one of the jets satisfies 
the tau candidate selection criteria. Another significant 
source of events with isolated muons and tau candidates 
is W + jets production, where the W boson decays to 
poj and one of the jets is misidentified as a tau candidate. 
The Z — > fi + [i~~ background is reduced by the require- 
ment that no other muon be found in the event, but a 
small number of events will be selected when one of the 
muons is not reconstructed. Small contributions are also 
expected from W — > tv and WW — > Ivlv, as well as 
ti production. Contributions from W Z and ZZ events 
yield less than one event each after the final selection 
criteria and are therefore neglected. All backgrounds, 
except the multijet background, are estimated using MC 
simulations. 

The multijet background is estimated using the data 
events that satisfy all requirements placed on the signal 
sample except that the muon and the tau candidate have 
the same sign charge. We call this the same-sign (SS) 
sample. To obtain the appropriate normalization for this 
background, a special data sample is selected, the "multi- 
jet sample," containing events that pass all requirements 
placed on the signal sample except the isolation criteria 
and the cut on the tau NN output. Instead of the iso- 
lation requirement used for the signal events, the events 
in the multijet sample have the sum of energies of all 
calorimeter cells inside the hollow isolation cone in the 
range 2.5 to 10 GeV. The sum of all non-muon tracks 
Pt within the track isolation cone is required to be in the 
same interval 2.5 — 10 GeV. To avoid contamination from 



7 



Z — > t + t~ signal events, an upper limit on the tau NN 
output is placed at 0.8. To increase the statistics of this 
sample, the muon p T is required to be at least 10 GeV 
instead of 15 GeV. The multijet sample is expected to 
be completely dominated by multijet processes, but may 
also include events in which a W decaying into a muon is 
produced in association with a jet. The W + jets contri- 
bution is reduced by requiring that the muon and the tau 
candidate are back to back (|(/> M — <p T \ > 2.5). A slight 
excess of opposite sign (OS) over SS events is observed 
in the multijet sample. No significant dependence of the 
OS/SS ratio as a function of pr or NN output is observed 
for the three types of tau candidates in the multijet sam- 
ple. Correction factors (f mj ) of 1.13 ± 0.03, 1.08 ± 0.01, 
and 1.06 ± 0.01 for tau types 1 to 3 are obtained, being 
used as discussed below to normalize the multijet back- 
ground in the final signal sample. 

The number of events in the SS sample is corrected 
for the contribution from Z —> Z — > t + t~, and 

W — > tv obtained from MC, multiplied by an additional 
correction factor which takes into account the difference 
between the charge misidentification rates in data and 
MC. Totals of 6 events for type 1, 16 events for type 
2, and 18 events for type 3 tau candidates from Z — > 
Z — ► t + t~, and W — ► tv are estimated to have a 
misidentified charge after all cuts and are subtracted from 
the number of SS events when the multijet background 
is calculated. The contribution from W — > fiv events is 
accounted for separately. 

A part of the W + jets background is already included 
in the SS sample. However, we expect a significant excess 
of OS events compared to the number of SS events due to 
the fact that a high percentage of W + 1 jet events comes 
from quark jets. The number of VF+jets events in data is 
estimated by selecting a sample that is expected to have 
a large contribution from W boson processes and low or 
negligible contributions from Z boson production. Such a 
VK+jets enriched sample can be obtained by requiring an 
isolated muon with px > 20 GeV, a tau candidate with 
0.3 < NN < 0.8, |</v -<t> T \< 2.7, and m* > 40 GeV. 
Mostly multijet and VF+jets events contribute to this 
sample. The excess of OS events compared to SS events 
is given for the multijet background by for tau type i. 
For the W + jets sample, similar factors (fyy) 01 2-39 ± 
1.01, 3.15 ± 1.17, and 1.6 ± 0.26 are estimated from data, 
in the sample with the cuts listed above, but requiring a 
tighter cut m* > 60 GeV. Using this, we can calculate the 
number of W + jets events in the W + jets enriched data 
sample by solving the following system of two equations 
for each tau type i: 



where N^ is the number of W - 



- jets events, N^- is the 



V, 



w + N^ = Nh s 



N, 



ss 



<w 



- N l 



f ■ - 1 

J mi 



N l ■ 
+ 1 mJ 



N, 



OS 



N 



ss 



(4) 



(5) 



number of multijet events and Nq S , Ng S are the num- 
bers of OS, respectively SS events in the W + jets en- 
riched data sample. The ratios between the number of 
W + jets events calculated in data by solving the above 
system of equations and the one expected from MC for 
each tau type are used as normalization factors for this 
background in the signal region. The uncertainty on N^ 
from data is taken as a systematic uncertainty. The esti- 
mated number of W + jets events in the signal sample, 
not including those in the SS sample, is 14 ± 5 events. 

Several distributions such as muon and tau candidate 
transverse momentum, pseudorapidity, and azimuthal 
angle, as well as ^7, m* , and visible mass are compared 
between the data and the predicted sum of backgrounds 
and Z — > t + t~ for the SM cross section and branching 
ratio. All these distributions show good agreement af- 
ter each of the preselection, NN selection, and anti-VF 
requirement stages. 

In Fig. [2] the visible mass distribution for events which 
pass the final selection requirements is shown separately 
for each of the tau types, while Fig. [3] shows the same 
distribution for the sum of all types. Good agreement 
is observed between the data and the sum of the back- 
ground SM processes and Z — ► t + t~ signal, normalized 
using the NNLO SM prediction 
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n. 



Table Q] shows the number of events expected for each 
tau type from each of the backgrounds, as well as from 
the Z t + t~ signal. It also shows the total numbers of 
expected background and signal events in comparison to 
the numbers of events observed in data for three levels 
of selection: preselection, preselection plus NN output 
requirement, and after all selection criteria are applied. 
Good agreement is observed between the predicted and 
observed numbers of events at each level of selection for 
all tau types. 

We estimate that approximately 1.2% of all Z — > t + t~ 
events have the wrong sign for either the muon or the 
tau candidate, therefore appearing as SS events. From 
the number of Z —* t + t~ events obtained by subtracting 
the estimated background from the number of events in 
the final sample, we calculate the number of Z — > t + t~ 
events reconstructed as SS to be 17. This number is 
added to the number of events in the OS sample when 
calculating the Z — * t + t~ cross section. 

Reconstruction of a second track close to a first recon- 
structed track is found to be more efficient in MC than 
in data. A correction factor of 0.97 ± 0.028 is applied to 
simulated events containing type 3 tau candidates. This 
factor is obtained by comparing the ratios of type 3 tau 
candidates with two and three tracks in data and MC 
and taking into account that there are twice as many SS 
as OS combinations when one of the three tracks is lost. 

Systematic uncertainties on the multijet and W + jets 
backgrounds are derived from the statistical uncertainties 
of the control samples used to estimate these backgrounds 
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FIG. 2: Visible mass distribution for (a) type 1 tau events, (b) type 2 tau events, and (c) type 3 tau events. The data are the 
points with error bars. The different components of the SM expectation are as given in the legend. The Z — > t + t~ signal is 
normalized to the theoretical expectation calculated at NNLO using MRST2004 PDFs [TH. fl^|. 





Type 1 


Type 2 


Type 3 


Process 


Preselection 


Preselection 


All cuts 


Preselection 


Preselection 


All cuts 


Preselection 


Preselection 


All cuts 






± NN > 0.9 






± NN > 0.9 






± NN > 0.95 






302 ±4 


230 ±4 


146 ±3 


1469 ±9 


1131 ± 8 


786 ±7 


693 ±6 


484 ±5 


358 ±5 


Zh* -> fi+fT 


58 ±2 


43 ±2 


6.1 ±0.6 


176 ±3 


108 ±3 


14.0 ±0.8 


184 ±3 


38±1 


8.9 ±0.7 


ww 


7.2 ±0.3 


6.1±0.3 


0.4 ±0.1 


79 ±1 


74 ±1 


6.9 ±0.3 


9.3 ±0.4 


6.1±0.5 


0.5 ±0.1 


ti 


2.7±0.3 


2.0 ±0.3 


0.2±0.1 


33 ±1 


28 ±1 


2.4±0.3 


29 ±1 


4.2 ±0.4 


0.5 ±0.1 


W -»TW 


10 ±2 


4±1 


1.5 ±0.8 


50±4 


14.1 ±2.2 


1.4 ±0.7 


168 ±7 


22.0 ±2.7 


3.7±1.2 


W -> \w 


127± 11 


42 ±5 


2.1 ±0.9 


470 ± 18 


116±9 


6.7 ±1.9 


1384 ±32 


202 ± 13 


14.1 ±2.7 


Multijet 


208 ±15 


46±8 


25±5 


584 ± 25 


123 ±12 


61±8 


2265 ± 47 


273 ± 18 


145 ±13 


Predicted 


715 ±18 


373 ±11 


181 ±7 


2861 ±32 


1594 ±18 


878 ±12 


4732 ± 59 


1029 ±23 


531 ±15 


Data 


720 


380 


170 


2836 


1546 


843 


4760 


981 


498 



TABLE I: Number of OS events expected for each tau type from the Z — > r + r~ signal as well as from each of the backgrounds, 
their sum and the number of OS events observed in data, for three levels of selection: preselection, preselection ± NN output 
> 0.9 (0.95 for type 3) and after all selection criteria are applied (preselection ± NN output > 0.9 or 0.95 ± m* < 20 GeV). 
The uncertainies are statistical. 



and from the systematic uncertainties on the correction 
factors used for their normalization. 

The systematic uncertainty related to the tau energy 
measurement is estimated by scaling the charged pion 
response used for data by the largest difference found 
between the response measured in data and the response 
obtained using gCALOR (6%) and recalculating the ac- 
ceptance applying all cuts. The value of this uncertainty 
is 1%. 

NN systematic uncertainties are calculated using sta- 
tistical ensembles of events in which each input variable 
is allowed to fluctuate within the difference observed be- 
tween the distributions of that particular variable in data 
and MC. The RMS of the ratio of the number of events 
passing a certain NN cut to the number of events in the 
ensembles, called the ensemble cut ratio, is taken as a 
measure of the NN uncertainty. The estimated uncer- 
tainties are 4.3% for type 1, 2.0% for type 2, and 3.8% 
for type 3 tau candidates, which results in a total uncer- 
tainty of 2.7%. 

The uncertainty due to the tau candidate track recon- 



struction efficiency is taken to be the same as the uncer- 
tainty on reconstructing muon tracks and is estimated 
using Z — ► events to be 1.4%. The uncertainty 

on the correction factor due to differences between data 
and MC in tracking efficiency for type 3 taus is added in 
quadrature to this value, resulting in a total uncertainty 
related to the tau candidate tracks of 1.6%. The system- 
atic uncertainties due to muon identification and muon 
track matching are determined to be 0.6% and 0.8%, re- 
spectively. The systematic uncertainty due to the charge 
misidentification is 1%. The uncertainty on the trigger 
efficiency is 2.7% and takes into account the bias related 
to the choice of the control sample, the variation due to 
possible background contamination, variations in time or 
due to changing luminosity, the choice of binning, and the 
choice of parameters for the efficiency, as well as the lim- 
ited statistics. The uncertainty on the total integrated 
luminosity is 6.1% [4], with an additional systematic un- 
certainty of 1% related to the influence on the luminosity 
of applying the data quality criteria used to reject events 
with coherent calorimeter noise. 



9 




150 200 250 

Visible Mass [GeV] 

FIG. 3: Visible mass distribution for all tau types. The 
Z — > t + t~ signal is normalized to the theoretical expecta- 
tion calculated at NNLO using MRST2004 PDFs. 



The PDF uncertainty of 2.0% is estimated using a NLO 
calculation [l9j and the CTEQ6.1 error sets. This uncer- 
tainty is obtained from the variation in acceptance when 
these error sets are used, added in quadrature with the 
difference in acceptance when using the MRST2004 er- 
ror sets at NLO and with the additional variation when 
going from NLO to NNLO with MRST2004. Table HJ 
summarizes all the systematic uncertainties. 



Source 


Value 


Tau energy scale 


1.0 % 


Tau identification 


2.7 % 


Tau track reconstruction 


1.6 % 


Multijet background 


1.6 % 


W — ► fiv background 


0.5 % 


Trigger 


2.7 % 


Muon track match 


0.8 % 


Muon identification 


0.6 % 


Muon momentum resolution 


0.4 % 


Charge misidentification 


1.0 % 


MC statistics 


0.6 % 


PDF 


2.0 % 


Total (except luminosity) 


5.2 % 


Luminosity 


6.2 % 



TABLE II: Systematic uncertainties on the u(pp — > Z/'y* 
X) ■ Br(Z/7* — > t + t~) measurement. 



The cross section times branching ratio for the process 
pp — > Z 77* + X — > t + t~ + X is given by the number 
of signal events divided by the product of the total ef- 
ficiency and the integrated luminosity. The number of 
signal events estimated from Table U with the correction 
for signal events reconstructed as SS, is 1227. Since Ta- 
ble Q] shows the estimated number of events in the Z/j* 
mass range 15 — 500 GeV, other corrections have to be 
made in order to compare the result of this analysis with 



theoretical cross sections. To limit the mass range to 
60 — 130 GeV, the number of events expected from the 
mass region 15 — 60 GeV (7 events) as well as from the 
130 — 500 GeV mass region (26 events) are subtracted 
from the number of signal events in data. The total effi- 
ciency for Z — ► t + t~ events in the 60 — 130 GeV mass 
region is 4.9 x 10 -3 , which also includes the trigger effi- 
ciency of 52.3%. Finally, a factor of 0.98 20] is applied 
to estimate the Z boson cross section as opposed to the 
Z/7* cross section for this mass region. 

Given the systematic uncertainties listed in Tabic HT1 
and an integrated luminosity of 1003 pb - , we estimate 
a{pp -> Z + X) ■ Br(Z -> r+r") = 240 ± 8 (stat) ± 
12 (sys) ± 15 (lum) pb, which is in good a g reement with 
the SM prediction of 251.9±^° 8 pb [H [jj that results 
from the NNLO calculation using the MRST2004 PDFs, 
as well as with the 241.6jl|;| pb [ll|, [HI value obtained 
at NNLO using the CTEQ6.1M PDF parametrization. 
This result is the most precise measurement of <j(pp — > 
Z + X) ■ Br(Z — > t + t~) to date, in good agreement with 
previous measurements of the Z boson cross section times 
branching ratio to leptons at y/s = 1.96 TeV @, H, [2l[. 
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